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Abstract 



A brief summary of results on homotheties in General Relativity is given, in- 
cluding general information about space-times admitting an r-parameter group 
of homothetic transformations for r > 2, as well as some specific results on 
perfect fluids. Attention is then focussed on inhomogeneous models, in partic- 
ular on those with a homothetic group H4 (acting multiply transitively) and 
H3. A classification of all possible Lie algebra structures along with (local) 
coordinate expressions for the metric and homothetic vectors is then provided 
(irrespectively of the matter content), and some new perfect fluid solutions are 
given and briefly discussed. 

"£5 1 Introduction 

cr 

This paper is devoted to the study of space-times admitting an intransitive group 
of homotheties, with a view towards those which can be interpreted as perfect fluid 
solutions of Einstein's field equations [fj|]. 
^ ■ A collection of important results regarding generic properties of space-times admit- 

ting homothetic transformations can be found in @-0 (and references cited therein), 
and in || where the case of multiply transitive action is thoroughly studied by Hall 
and Steele. 

The study of this subject began with the pioneering paper by Cahill and Taub 
||, followed by the works of Eardley |T0| , |TT|j . From then on, homotheties have 



been studied in connection with a wealth of situations of physical interest in classical 



general relativity as well as in cosmology, see fl2| , [L3[ |14|| for interesting reviews on 
homothetic solutions. 

The paper is organized as follows: section 2 contains a brief summary of results 
on groups of homotheties and space-times admitting them, the implications that 
they have on perfect fluids and the physical quantities characterizing them (density, 
pressure, velocity,...), and we also summarize all the general information about space- 
times admitting an r-parameter group of homothetic transformations for r > 2. This 
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includes: dimension of the nomothetic and isometric algebras as well as that of the 
orbits they act on respectively, together with its nature (spacelike, timelike or null), 
allowed Petrov and Segre types of the Weyl and Ricci tensors, and whether perfect 
fluid solutions exist or not. Most of the contents of this section is a review of dispersed 
results in the literature, but we considered useful to gather them all in a single table. 
Some of the results are, as far as we are aware of, new (see especially those cases 
where null orbits occur). 

In sections 3 and 4 a classification of all possible Lie algebra structures along with 
(local) coordinate expressions for the metric and homothetic vectors, are given for the 
cases H4 (acting multiply transitively on three-dimensional orbits, in keeping with our 
assumption of intransitive action) and H3 respectively. These classes of space-times 
can be understood as generalizations to the Kantowski-Sachs and Bianchi models 
respectively. The characterizations provided are independent of the field equations, 
and therefore they may have applications other than those considered here (perfect 
fluids) . 

In particular, section 3 contains some review material on the case r = 4, together 
with some new results, all of them presented in a unified manner, extending the work 
of Wu [|1^] , Cahill and Taub and Shikin ]16| . The general perfect fluid solution is 
then given in certain, well-defined and invariantly characterized subcases. Whenever 
this is not possible, a few selected examples are presented. It is assumed that the 
matter satisfies the weak and dominant energy conditions, and expressions for the 
kinematical quantities (acceleration, expansion, deceleration parameter, shear and 
vorticity) are provided for each case. 

Finally, section 4 contains some new solutions for the case r = 3 and appropriate 
references to related work on this issue. We summarize the results concerning the 
topology of the Killing orbits and the Bianchi classification of the homothetic algebras. 
We distinguish the cases where the Killing subalgebra is Abelian from that where it is 
non-Abelian. Attention is then restricted to the orthogonally transitive case, giving 
for each possible Lie algebra structure the coordinate forms of the proper homothetic 
vector field and the metric. In the Abelian case we distinguish three different classes of 
such models, depending on the orientation of the fluid flow relative to the homothetic 
orbits. The case in which we are more interested is the so-called "tilted" where new 
solutions are found. Finally, we provide explicit forms for the homothetic vector field 
and the metric in the case of a (maximal) non-Abelian G2, and, although no perfect 
fluid solutions have been found, we briefly discuss some properties. 



2 Basic facts about homotheties 
2.1 Definition and properties 

Throughout this paper (M, g) will denote a space-time: M then being a Hausdorff, 
simply connected, four-dimensional manifold, and g a Lorentz metric of signature 
(-,+,+,+). All the structures will be assumed smooth. 

A global vector field X on M is called homothetic if either one of the following 
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equivalent conditions holds on a local chart 



£x9ab = 2ng ab , X a . b = ng ab + F ab , (1) 

where n is a constant on M , £ stands for the Lie derivative operator, a semi-colon 
denotes a covariant derivative with respect to the metric connection, and F ab = —F ba 
is the so-called homothetic bivector. If n 7^ 0, X is called proper homothetic and if 
n = 0, X is a Killing vector (KV) on M. For a geometrical interpretation of ([l]) we 
refer the reader to |2|. |6| . 

A necessary condition that X be homothetic is 

A ;6c — -K bcd^ , (/J 

where i? a b c <i are the components of the Riemann tensor in the above chart; thus, 
a homothetic vector field (HVF) is a particular case of affine collineation |7j and 
therefore it will satisfy 

£-xR a bcd = C-xRab = £xC a bcd = , (3) 

where R ab (= R c a cb) an d C a ft c d stand, respectively for the components of the Ricci 
and the conformal Weyl tensor. 

The set of all HVFs on M forms a finite dimensional Lie algebra under the usual 
bracket operation and will be referred to as the homothetic algebra, 7i r , r being its 
dimension. The set of all KVs on M also forms a finite dimensional Lie algebra, the 
Lie algebra of isometries, which will be denoted as Q s (s being its dimension), and 
one has that Q s C 7i r (i.e., Q s is a subalgebra of TC r ). Furthermore, it is immediate 
to see by direct computation that the Lie bracket of an HVF with a KV is always a 
KV and that, given any two proper HVFs, there always exists a linear combination 
of them which is a KV. From these considerations it immediately follows that the 
highest possible dimension of 7i T in a four- dimensional manifold is r = 11. 

If r 7^ s then s = r—1 necessarily, and one may choose a basis {Xx, • • • , X r -i, X} = 
{Xa}a=i---t for Tir-, m such a way that X is proper homothetic and X\, • • ■ , X r ^\ are 
Killing vector fields spanning Q r ~\- If these vector fields in the basis of 7i r are all 
complete vector fields, then each member of 7i r is complete and Palais' theorem 
Pj, [17], |H| guarantees the existence of an r-dimensional Lie group of homothetic 
transformations of M (H r ) in a well-known way; otherwise, it gives rise to a local 
group of local homothetic transformations of M and, although the usual concepts of 
isotropy and orbits still hold, a little more care is required H. 



The following result || IS] will be useful: 

The orbits associated with TC r and Q r -i can only coincide if either they are four- 
dimensional or three-dimensional and null. (This result still holds if 7i r is replaced 
by the conformal Lie algebra C r and does not depend on the maximality of 7i r or C r ). 

The set of zeroes of a proper HVF, i.e., {p G M : X(p) =0} (fixed points of the 
homothety), either consists of topologically isolated points, or else is part of a null 
geodesic. The latter case corresponds to the well-known (conformally flat or Petrov 
type N) plane waves j2|, [20f . 
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At any zero of a proper HVF on M all Ricci and Weyl eigenvalues must necessarily 
vanish and thus the Ricci tensor is either zero or has Segre type {(2, 11)} or {(3, 1)} 
(both with zero eigenvalue), whereas the Weyl tensor is of the Petrov type O, N or 
III |2| (see also |2l| for vacuum space-times). 



2.2 Perfect fluids 

The energy-momentum tensor of a perfect fluid is given by 

Tab = (a* + P)u a u b + V9ab , (4) 

where /i and p are, respectively, the energy density and the pressure as measured by 
an observer comoving with the fluid, and u a (u a u a = — 1) is the four- velocity of the 
fluid. If X is an HVF then, from Einstein's Field Equations (EFE) it follows that 

C-xT a b = , (5) 



and this implies in turn |T0| 



C x u a = nu a , C x p = -2np , C x ^ = -2n// . (6) 

Thus, the Lie derivatives of u a , p and fi with respect to a KV vanish identically. 

If a barotropic equation of state exists, p = p(fi), and the space-time admits a 
proper HVF X then [fZ| 

p=( 7 -l)/i, (7) 

where 7 is a constant (0 < 7 < 2 in order to comply with the weak and dominant 
energy conditions). Of particular interest are the values 7 = 1 (pressure- free matter, 
"dust") and 7 = 4/3 (radiation fluid). In addition, the value 7 = 2 (stiff- matter) 
has been considered in connection with the early Universe. Furthermore, values of 7 
satisfying < 7 < 2/3, while physically unrealistic as regards a classical fluid, are 
of interest in connection with inflationary models of the Universe. In particular, the 
value 7 = 0, for which the fluid can be interpreted as a positive cosmological constant, 
corresponds to exponential inflation, while the values < 7 < 2/3 correspond to 



power law inflation in FRW models [ 23"| , but it is customary to restrict 7 to the range 
1 < 7 < 2. 

If the proper HVF X and the four- velocity u are mutually orthogonal (i.e., u a X a = 
0) and a barotropic equation of state is assumed, it follows that 7 = 2, i.e., p = ji 



stiff-matter |10|| , on the other hand, if X a = au a the fluid is then shear-free. Further 



information on this topic can be found in pi, [25|, |26| 



2.3 The "dimensional count-down" 

In this subsection, the maximal Lie algebra of global HVF on M will be denoted as 
7i r (r being its dimension), and it will be assumed that at least one member of it is 
proper homothetic. 

The case of multiply transitive action is thoroughly studied in ||, and we shall 
refer the reader there for details; nevertheless, and for the sake of completeness, we 
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summarize in the following table the results given there, which follow invariably from 
considerations on the associated Killing subalgebra and the fixed point structure of the 
proper HVF. Futhermore, we have added a few other results, also in the literature or 
following straightforwardly from those, so as to complete the study down to dimension 
4. 
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Table 2.1 



The first entry in the table gives the dimension of the group of homotheties, the 
second and third entries stand for the nature and dimension of the homothetic and 
Killing orbits respectively (e.g.: N 2 , T 2 and S 2 denote Null, Timelike and Spacelike 
two-dimensional orbits respectively, 3 stands for three-dimensional orbits of either 
nature, timelike, spacelike or null), the fourth and fifth entries give the Petrov and 
Segre type(s) of the associated Weyl and Ricci tensors (in the latter case it is to 
be understood that all possible degeneracies of the given types, can in principle oc- 
cur, including vacuum when possible). Finally, the last two entries give respectively 
the possible interpretation whenever it is in some sense unique, and the existence 
or non-existence of perfect fluid solutions for that particular case, along with some 
supplementary information; thus FRW stands for Friedmann-Robertson- Walker, LRS 
for Locally Rotationally Symmetric, and Bianchi refers to that family of perfect fluid 
solutions. The cases that cannot arise are labeled as "Not Possible", and wherever 
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no information is given on the Petrov and Segre types, it is to be understood that 
all types are possible in principle. The Segre type of the Ricci tensor of the case de- 
scribed in the last row, is unrestricted except in that it must necessarily have two equal 
(spacelike) eigenvalues; perfect fluid solutions of these characteristics constitute spe- 
cial cases of spherically, plane or hyperbolically symmetric perfect fluid space-times. 
For further information on LRS space-times, see [27, 28]; for the case r 
and null three-dimensional Killing orbits, see 

|9 



29 



neous Bianchi models, see 



35, 36 



4 transitive 
Regarding spatially homoge- 
and for the last three cases occurring in 
and II. 



the table, see respectively [31j, |32 

The case r = 3 has an associated Killing subalgebra Q2 and the respective dimen- 
sions of their orbits are 3 and 2 



(see for instance |33], [M], [42|] and references cited 
therein). When the Killing subalgebra has null orbits, the metric is of Kundt's class 
and perfect fluids are excluded. If the Killing orbits are timelike, the solutions 
can then be interpreted as special cases of axisymmetric stationary space-times (pro- 
vided that regularity conditions hold on the axis JI|, |3~9fl), and if they are spacelike 
as special cases of inhomogeneous cosmological solutions or cylindrically symmetric 
space-times. In both cases, perfect fluid solutions have been found. 



3 The H4 case 

The associated isometric group for perfect fluid space-times acts necessarily on T3, 
S3 or S2 orbits (see Table 2.1). In the intransitive case, the G3 must act multiply 
transitively on two-dimensional surfaces of maximal symmetry S2, which are then of 



constant (positive, zero or negative) curvature and admit orthogonal surfaces [HJ . 

Possibly no problem in this context has been more exhaustively studied than that 
of spherically symmetric homothetic space-times, beginning with the seminal paper 



of Cahill and Taub || and continuing with recent papers by Ori and Piran |45fl , Carr 



and Yahil f4"6| , Henriksen and Patel j|7||, and Foglizzo and Henriksen |58| among 
others (see references therein). Homothetic space-times with plane symmetry are 
also considered by Shikin |T6f . 

What we attempt in this section, rather than presenting a survey of the models 
existing in the literature, is to study in a unified manner all possible cases, i.e., 
homothetic orbits of either nature (timelike, spacelike or null at every point peM 
and the more general case in which their nature varies from point to point) as well 
as the different possibilities for the curvature k of the isometry orbits (spherical and 
plane symmetry as well as the k = — 1 case); thus, extending previous works by Wu 



T5| , where only spacelike homothetic orbits are considered (i.e., type B and some 
type C solutions in our classification below), and by Cahill and Taub [[|, and Shikin 
IB] where only spherical and plane symmetry are considered respectively. 
As it is well known, the space-time metric can be written as [p]]: 

ds 2 = -A 2 (r,t)dt 2 + B 2 (r,t)dr 2 + F 2 (r,t)(dy 2 + Y 2 (y,k)dz 2 ) , (8) 
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{sin y k — +1 
y k = (9) 
sinh ?/ k = —1 . 

The Killing vectors being ^ = sin^c^ + ^ cos zd z , £2 = coszd y — ^ sin 2:<9 2 , and 
6 = d z , where the dash denotes a derivative with respect to y, and they satisfy the 
following commutation relations 

[6,6] = ^6, [6,6] =6, [6,6] =6- (10) 

Assuming the existence of a proper HVF, X, and since its commutator with a KV 
must be a KV, the Jacobi identities imply the following structures for 7^ 4 

(/) [x,6] = o, j = 1,2,3, fc = o,±i, (11) 

X = X t (t,r)d t + X r (t,r)d r , (12) 
(J/) [X, £1] = ^ , [X, £ 2 ] = 6 , [X, 6] = , fc = , (13) 
X = X% r)d t + X r (t, r)d r - 1/5, . (14) 

Requiring the coordinate system to be a comoving one (i.e., u a = —A{t,r)5 t a ) the 
HVF takes the form 

X = X\t)d t + X r (r)d r + Xy(y)d y = X + Xy(y)d y , (15) 

and the following possibilities then arise: 

{A) X = d t , (16) 

(B) X = d r , (17) 

(C) X = d t + d r . (18) 

The form (A) corresponds to u being tangent to the timelike homothetic orbits, T%, 
(B) corresponds to the case of spacelike homothetic orbits, S3, orthogonal to the fluid 
flow, and (C) is the most general (tilted) case, including also the possibility of having 
null homothetic orbits, A3, when the functions A(t, r) and B(t,r) in @ equal each 
other. 

The homothetic equation ([!]) specialized to the metric (H) yields then the following 
possibilities: 



Case (A) 



A 


k 


X 


A 2 (t,r) 


B 2 (t,r) 


F 2 (t,r) 


I 


-1,0,1 


d t 


e 2t H 2 (r) 


e 2t H 2 {r) 


e 2t f 2 {r) 


II 





d t - ydy 


e 2nt H 2 {r) 


e 2nt H 2 {r) 


e 2[n+l)t 



Table 3.1 
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where t has been scaled in A I so as to make n = 1 in ([!]). Solving now the field 
equations for a perfect fluid in each one of the two above cases, we find that the only 
possible solution for the family AI has an equation of state \x + 3p = and it admits 
(at least) one further Killing vector, thus being a particular case of TC5 and we shall 
not study it here. 

For All, two families arise which depend on the value of n: 



All 


H(r) 


f(r) 




7 


n G (— 00, — 3)U 
(-2,-l)U(0,+oo) 


a(cosh(ar))" 


(f) n+1 


P p -2nt rr2(a-l) 

a* 5 


2n 


2 + 3n 


n e (-3,-2) 


a(sinh(ar))~~ 


(¥) n+1 


P p -2ntTT2(a-l) 


2n 


2 + 3n 



where 



Table 3.2 

2(n +1) (n + l)(n + 3)(3n + 2) 



_ -r Ljyu -r -r ^ 

n + 2 + 2) 

and a is a constant. The vorticity is zero and the volume expansion 9, deceleration pa- 
rameter q = — 1 — 39/9 2 , acceleration ii a , and non- vanishing shear tensor components 
a a b can be given as: 
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2 + 3n 

< 

H 



2e nt H 



-nt 



Pi, 



q 



2 + 3n 

e (2+n)i p 



U 



e (2+n)t y 2 J2 



(20) 



[21\ 



Ml 3H 
Notice that, depending on the value of n , the solution contracts and decelerates or 
it expands and inflates. 

With regard to the dimensionless scalars, the density parameter Q = 3n/9 2 , the 
dimensionless acceleration W = ii/9 and the shear parameter S = 3a 2 /9 2 (see [^1|] 
for further details), one has for the first family 



n 



3(n + l)(n + 3) 



1 



(2 + 3n)n 2 (n + 2)cosh 2 (ar) 



W 



1 



2 + 3n 



tanh(ar) 



(221 



the models then being accelerated dominated at large distances; whereas for the 
second family 



n 



3(n + l)(n + 3) 



(2 + 3n)n 2 (n + 2)sinh 2 (ar) ' 



W 



1 



2 + 3n 



coth(ar) 



(23) 



thus being asymptotically spatially homogeneous. In both cases lim^oo Q = 0, (vac- 
uum dominated models) and E = 1/(2 + 3n) 2 , as this quantity is non-vanishing for 
all possible values of n, the models have no isotropic limit. 



Case (B) 
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B 


k 


X 


A 2 (t,r) 


B 2 (t,r) 


F 2 (t,r) 


I 


-1,0,1 


d r 


e 2r H 2 (t) 


e 2r H 2 (t) 


e 2r f(t) 


II 





d r - ydy 


e 2nr H 2 {t) 


e 2nr H 2 (t) 





Table 3.3 



where r in case BI has been re-scaled so as to have n — 1. 

Solving now the field equations for a perfect fluid source, one has: 



BI 








Domain 


k = ±1 
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e 2 ' + § 


4e 2r f 4 


t : e 2 * + | > 


fc = 0,±1 
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a sinh 2t + § 


A; 2 + 4a 2 
4e 2r .f 4 


* : a sinh 2t + | > 


k = ±1 


1 


a cosh 2t + | 


A; 2 - 4a 2 
4e 2r / 4 


t : a cosh 2t + | > 



Table 3.4 

where a is an arbitrary constant, which must be different from zero to prevent the 
occurrence of further Killing vectors, k ^ in both the first and last cases since 
otherwise one would have a vacuum solution in the former case, and negative energy 
density in the latter. In all three cases, the fluid is irrotational and has an stiff 
equation of state, i.e., 7 = 2. For these solutions one has 



9 = 



Q = 



2e r f 



3 // 



2 



u 



— d r 



W 



f_ 
2/ 



a yy 



e r ff 



e r E 2 // 



(24) 
(25) 



3/ ' ~ yy 3 ' 3 

a dot indicating a derivative with respect to t. In all cases the shear parameter is 
£ = 1/4, thus there is no isotropic limit. For all solutions in BI lim^oo W = 1/2, 
thus corresponding to accelerated dominated models, hence with no FRW limit. 
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r\t) 




n < -3 
n > -1 


a[S'] 2 ("+ 1 ) 
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(n + l)(n + 3) e _ 2nr[s] -^ 


-3 < n< -1 


a[C] 2 ("+D 
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(n + l)(n + 3) e _ 2nr[cr ^ 
a 2 n 2 L 1 



Table 3.5 



where 



a = constant , S = sinh(2(n + l)t) , C = cosh(2(n + l)t) . (26) 
In both cases the fluid is irrotational, its equation of state being p = fi, and one has: 



u = nd r , 6 = e~ nr H- x + 2 f -. 
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n 



(2n + 3) 2 



(27) 



As for the respective dimensionless scalars 

3(n + l)(ra + 3) 1 



n = 

Q = 

for the first case, and 

n = 
q = 



w 



n 



n 2 (2n + 3) 2 C 2 ' " 2n + 3 C 

-(An 2 + 6n + 1) C 2 + 6(w + \){2n + 3) 
(2n + 3) 2 C 2 

3(n+ l)(n + 3) 1 



n C 



n 2 (2n + 3) 2 S 2 ' " 2n + 3 S 
(4n 2 + 6n + 1) S 2 + 6(n + l)(2n + 3) 



(28) 
(29) 

(30) 
(31) 



(2n + 3) 2 S 2 

for the second case. So, the solutions behave in different ways depending on the value 
of the parameter n and the range of values of t considered. 

Case (C) 



c 
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A 2 (t,r) 


B 2 (t,r) 


F 2 (t,r) 
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-1,0,1 


d t + d r 


e t+r H 2 (t - r) 


e t+r L 2 (t-r) 


e t+r f 2 (t - r) 


II 





d t + d r - yd y 


e n(t+r ^H 2 (t-r) 


e n ^L 2 (t-r) 


e (n+l)(t+r) p( t _ r ) 



Table 3.6 

where t and r in CI have been re-scaled so as to make n — 1. 

Unfortunately, no general solutions to the field equations for a perfect fluid can be 
given in these cases, and as it turns out when trying to solve them in particular cases, 
most of the solutions thus found hold only on some open domains of the manifold of 
the form t — r > constant, nevertheless we next give two solutions that are valid over 
the whole space-time manifold, both of them corresponding to the case CI I: 



(CII, n = -3) 
Then one has 



/ = !, H 2 = ae^- 4 ^ , L 2 = ae^ . 
ll=v = e 3 (^) Q; - 1 [(4 - /3) e ( 4 -W*-'-) + (3e-^\ , 



(32) 
(33) 



where a and f3 are constants which can easily be chosen so as to make /i > all over 
M. The dimensionless scalars are 
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(f3 - 7) 2 
W = 



4- + (3e- i{t - r) 

1 -P -2(t-r) 

P-7 



, e 



P-7 



W - 7) 2 ' 



(34) 
(35) 
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(CII, n = 1) 



./ 



and then 



H 2 = aL 2 



L 



exp 



(1 _ a){t - r ) - W ~ <*) e -W-r)/(l-a) 



(fi-p)e t+r H 2 
(ft + p)e t+r H 2 



4(1 -«) , 
2(1 -a) 



[1 - a) + pe 



-2(t-r)/(l-a) 



(36) 



(37) 
(38) 



where again a and (3 are constants. The energy conditions restrict a to values < 
a < 1, and if we demand that the solution be valid over the whole manifold, then (3 
must be positive. Notice that in this case, for (3 ^ 0, there is no equation of state of 
the form p = p(fi). For (3 = the solution is an special case of Ti^. The dimensionless 
scalars are 



n 



3(1 -a) 



a 



w 



-1 + 



3-| + f e -I^(*-r> 



6 



^ - 



-1 



6 



1 2 



6 — a + (3e 1 ~ a 



(t-r) 



(39) 



(40) 



6-a + pe~— {t ~ r) 

In both cases there is no FRW limit. 

As a final remark to this section, notice that expressions appearing in tables 3.1 
3.3 and 3.6 are completely general, i.e., valid regardless the material content. 



4 The Hi case 



531], correct some errors and present new 



In this section we extend a previous work 
solutions. 

The existence of a 3-parameter homothetic group H 3 , implies that of a G2 C H 3 of 
isometries being the dimensions of their orbits three and two respectively (see section 



2.1). When the Killing subalgebra has null orbits, the metric is of Kundt's class [[E 
and perfect fluids are excluded [ 
Killing orbits are non-null. 



We shall therefore assume in the sequel that the 
A classification of all such space-times in terms of the Bianchi type of the ho- 



mothetic algebra can be found in [Q. See also the references cited therein for a 
(partial) account of papers on this issue. We can summarize the results concerning 
the topology of the Killing orbits and the Bianchi type of H 3 in the following table: 



G 2 - type 


G 2 — orbits 


i?3 — Bianchi type 


G 2 I 


S 1 x IR 


/, II, III 




I, II, III, IV, V, VI, VII 


G 2 II 


IR 7 


III 



Table 4.1 
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As the above table shows, two different topologies are possible in the Abelian G 2 
case (G 2 I) []33] namely, V 2 diffeomorphic to R 2 or to S 1 x 1R; and it follows in the 
latter case that the only possible Bianchi types for 7i 3 , irrespectively of the assumed 
matter content, are J, II and /// (this holds also if the is replaced by a conformal 
algebra C3, see ]39|]); and for the case V 2 = H 2 , the seven soluble Bianchi types can 
occur. For the non- Abelian case the only possible homothetic algebra is of the Bianchi 
type III, and its orbits are diffeomorphic to IR 2 . 



4.1 Case Gi Abelian 



In this subsection we shall restrict ourselves to the Abelian case with spacelike iso- 
metric orbits diffeomorphic to IR 2 , giving appropriate "translation rules" for the other 
possibilities. Furthermore, we shall assume that the Killing orbits admit orthogonal 
two-surfaces (i.e., orthogonally transitive G 2 metrics). Cosmological models admit- 
ting an Abelian G 2 on spacelike orbits have been studied by Ruiz and Senovilla [|50| 



and Van den Berg and Skea |52]| among others. The non-orthogonally transitive case, 
Wainwright's classes A(i) and A(ii), have been studied in []53| and []53J respectively. 

Adapting two coordinates to two commuting KVs, say £ = d x and r\ = d y , and 
choosing two other coordinates, t and z, on the surfaces orthogonal to the isometry 
orbits; it follows that the line element can be written in the form (see for instance 



H) 



ds 2 



-Adt 2 + Bdz 2 + R 



F(dx + Wdyf + F~ x dy' 



(41) 



where A , B , R , F , and W are all functions of t and z alone. 

All the other cases (i.e., timelike Killing orbits and Killing orbits diffeomorphic 
to S 1 x 1R of either nature, spacelike or timelike) can be formally obtained from the 
above by means of the following substitutions: 



T 2 


^ 1R 2 


V 2 S 1 x 1R 


d x 


h-> d t 




(t,x) H 


-* i(—x, t) 


Regularity condition 


W 


^ iW 


on the axis 



Table 4.2 

where <p is the angular coordinate (with the standard periodicity 2n). Regarding 
solutions with an Abelian G 2 acting on timelike orbits, (including the astrophysically 
relevant stationary and axisymmetric models, which have been studied for many 



years), it is worth mentioning that they have attracted renewed attention, see [[37 



It is easy to see from the commutation relations of the proper HVF, X, with £ 
and rj, and the homothetic equation specialized to the components g tx , 9ty, 9zx and 



g zy of the metric fl4TD that, X must take the form 

X*(t, z)d t + X z (t, z)d z + X x (x, y)d x + X y (x, y)d y 
X + X x (x,y)d x + X y (x,y)d y , 



X 



(42) 
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where X x (x, y) and X y (x,y) are linear functions of their arguments that yield for 
every different Bianchi type the following forms: 



Tvoe 


X x (x v) 


X y (x v) 


I 








II 


y 





III 


X 





IV 


x + y 


y 


V 


X 


y 


VI 


X 


qy 


VII 


-y 


x + qy 



Table 4.3 



As we are interested in perfect fluid solutions for the metric (f4lD, it is always 
possible to perform a change of coordinates in the t, z plane so as to bring the four- 
velocity of the fluid to a comoving form, preserving the diagonal form of the metric 
)5[. As a consequence, the fluid flow velocity u can be written as 



u 



l_d_ 
Adt 



or equivalently 



A o,o,o) , 



(43) 



the Einstein's field equations taking then a much simpler form. 

In this comoving coordinate chart, and taking into account the first equation in 
(|jD, it is easy to see that the part of the homothetic vector field orthogonal to the 
Killing orbits, X, is 

X = X t (t)d t + X z (z)d z . (44) 



We can now use the remaining coordinate freedom in the t, z plane [t — > m(t) , z 
n(z)} to bring X to either of the following three forms 



{ii) 
(in) 



X 
X 
X 



t ) 



d t + d r 



(45) 
(46) 
(47) 



Thus, three classes of perfect fluid solutions arise, depending on the orientation 
of the fluid flow u relative to the homothetic orbits: 

(i) The fluid flow is tangent to the homothetic orbits, and they are 
then timelike. 

(ii) The fluid flow is orthogonal to the homothetic orbits, and therefore 
they are spacelike. 

(Hi) "Tilted" fluid flow, i.e., u is neither tangential to nor orthogonal 
to the homothetic orbits, which are then not constrained -a priori- to 
being timelike or spacelike, and so that their nature may vary from point 
to point over the space-time. 
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4.1.1 Fluid flow tangent to the homothetic orbits 



This case, assuming the existence of two hypersurface orthogonal KVs (i.e., diago- 
nal metric) has been thoroughly studied by Wainwright, Hewitt, and collaborators 
in an interesting series of articles flUj |42l [431 , where the properties of these mod- 
els are analyzed using the qualitative theory of plane autonomous systems, showing 
that (first-class) self-similar solutions within this family can represent the asymptotic 
states at later times of more general inhomogeneous G2 models. Uggla |33 found four 
explicit solutions of this type. We shall not give here any explicit solution belonging 
to this class (see the above references and those cited therein), but rather provide the 
general form of the metric functions (including the non-diagonal cases) and briefly 
discuss the generic behaviour of the kinematical quantities associated with the fluid 
(acceleration, deceleration parameter, shear,...). 

The metric (in comoving coordinates) for the case of four-velocity tangent to the 
H 3 orbits (i.e., form (i) of the proper homothetic vector field) takes the form: 

ds 2 = e 2nt {-A(z)dt 2 + B{z)dz 2 + e at R{z)[F{t, z)(dx + W(t, z)dy) 2 + F~\t, z)dy 2 }} 

(48) 

where n is the homothetic constant. By rescaling the coordinate z one can set 



A(z) 



B(z) 



(49) 



Then, the functional form of the metric functions can be worked out for each Bianchi 
type. Thus, one has: 

(50) 
(51) 
(52) 
(53) 
(54) 

e {1 - q)t w(z) , (55) 
(56) 



CO 


a = 


. 


F = 


f(z) , W = w(z) , 


(II) 


a = 


0, 


F = 


f(z), W = w(z)-t, 


(III) 


a = 


-1 , 


F 


= e- l f(z) , W = e*w(z) , 


(IV) 


a = 


-2 , 


F 


= f(z), W = w(z)-t, 


(V) 


a = 


-2 , 


F 


= f(z), W = w(z), 


(VI) 


a = 


-(1 + 9) 


, F = e-^ t f(z), W = 


(VII) 


a = 


-s , 







w 



? + 



^/l + c(z) 2 + g(z) 2 + c{z) cos(y / 4 - q 2 t) + g(z) sm(^4 - q 2 t) 

^P^[-g(z) cos(y/4^t) + cjz) sin(y/4^fo)] 
1 + c(z) 2 + g(z) 2 + c(z) cos(V4 - qH) + g(z) sm(^A - q'H) 



Notice that the form of these functions again holds for any energy-momentum tensor, 
since no use has been made of the field equations in deducing them. 

The cases studied by Wainwright and collaborators correspond to types /, 
V, and VI, since these are the only ones in which the function W can be set equal to 
zero. For type VII, W — | implies the existence of a further Killing vector tangent 
to the Killing orbits and the metric would then admit a multiply transitive group 
of homotheties. Notice also, that the diagonal cases are separable in the variables t, 
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z, thus being special cases of the solutions studied by Ruiz and Senovilla J5(J or those 
of Agnew and Goode |5T] for 7 = 2. 

Specializing equation (^) to the matter variables fi and p, we obtain 



-2nt 



p 



~2nt 



p{z) 



(57) 



and by computing the kinematical quantities associated to the fluid velocity vector 
for the metric 



e 



3n + a 



ollt 



A 



q 



—a 



3n + a 



2~X 



dz 



(5* 



where a dash denotes a derivative with respect to z and q here denotes the deceleration 
parameter. The non-vanishing components of the shear tensor are 



a, 



—e 



Tit 



a 



,«* RF ( 3 E_ 



6VA 



'fit 



&xy ^ 



3 ' 

RF ( / ' • N 

3W- + Wa + 3W 
6y/A \ F , 



+ a 



D nt R 



a 



tin 



^ , -3- + a + 3W 2 FF + W 2 F 2 a + QF 2 WW 



(59) 

(60) 
(61) 



thus, the shear scalar is 



a' 



3(f) 2 + a 2 + 3F 2 W 2 



12Ae 



2nt 



(62) 



If one assumes that the fluid has an equation of state of the form (^), from the 
contracted Bianchi identities it follows 



7(/i + fx9)u a + 'jfxu a + (7 - l)n,b9 



ba 







where (1 = n^u*. Contracting the above expression with u a , one gets 



(63) 



(64) 



Assuming 6^0 (the case 6 = 0, although mathematically possible, is not physically 
interesting since it would correspond to a non-expanding universe, thus contradicting 
observations) and substituting (|57| ) and (|58D in equation (|64]), it follows 



7 



2n 



3n + a 

Specializing the quantities 7 and q to each Bianchi type, one gets 



(65) 



Type 


/ 


i7 


III 


IV 


V 


VI 


VII 


7 


2 
3 


2 
3 


2n 


2n 


2n 


2n 


2n 


3n-l 


3n-2 


3n-2 


3n-(l + g) 


3n — q 


q 








1 


2 


2 


1 + q 


Q 


3n-l 


3n-2 


3n-2 


3n-(l + q) 


3n — q 
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Table 4.4 



Notice that the only shear-free solution is the type I one. For types II to VI 
there is no limit where S = 3cr 2 /6 2 becomes zero; and for type VII, X — > for some 
spatial limit if and only if c(z) — > 0, g(z) — > and g = (i.e., F — > 1 and — > 0). 
Thus, only types I and V// can have solutions with a FRW limit, but in both cases 
fi + 3p = and q = 0, so they are not very relevant from a physical point of view. 



4.1.2 Fluid flow orthogonal to the homothetic orbits 



The case of spatially homothetic orbits was thoroughly studied by Eardley |T(J where 
a classification scheme of these models was given and their dynamical properties were 
studied. Luminet 



56 



constructed a convenient basis of 1-forms and gave its explicit 
form in terms of a standard coordinate basis {dx a } as well as the expression of the 
homothetic vector in the dual basis {d/dx a }. He also proved a theorem showing that 
perfect fluid models of a certain class were incomplete in the sense of Hawking and 
Ellis p . 

For the sake of completeness, we just mention that the form of the metric, corre- 
sponding to a homothetic vector field of the form (ii) , can be formally obtained from 
fl48|) to (|56|) by simply reversing the roles of the coordinates t and z. 

The expressions of the acceleration, expansion, deceleration parameter and shear 
scalar are given by 




u = ndz , 9 = -= \ — + 2- \ , (66) 



*(4) 2 + 4i + *(S) 2 -3|-6 



R 



(4+2$ 

3(f) 2 + 3F 2 (W/) 5 



(67) 



° = 12A^ • (68) 

In this case, the homothetic vector field, X, and the four- velocity u are mutually 
orthogonal, thus if a barotropic equation of state is assumed, then necessarily p = fi, 
i.e., stiff- matter and by simple inspection of the field equations this is seen to be 
equivalent to: 



| = (2n + a) 2 



(69) 



4.1.3 "Tilted" fluid flow 

Finally, the form (Hi) for the proper homothetic vector field is precisely the case we 
are currently interested in, namely, u not tangent nor orthogonal to the homothetic 
orbits. 
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Specializing now the nomothetic equation to the metric (41) we obtain 



ds 2 = e n{t+z) {-A(t- z)dt 2 + B{t- z)dz 2 

+e at ^R(t - z) [F{dx + Wdyf + F~ l dy 2 ] } , (70) 

where again n is the homothetic constant. The parameter a and the functional form of 
the metric functions for each Bianchi type, are those given by (|50~D - ([56|) after effecting 
the following substitutions 

t,— >^±£, z ^t-z. (71) 
The kinematical quantities being 




3n + a + f + 2f 



(72) 



2 _ (f-f + f) 2 + 3(») 2 + 3FW _ 

I 2e n{t +Z)A > V 6 ) 

where a dash denotes here a derivative with respect to (t — z). A careful study shows 
that there are no shear-free solutions in this case admitting a maximal group H 3 of 
homotheties: shear-free solutions are not possible for types II and IV; for the Bianchi 
types /, V and VII, the shear-free condition implies that the functions F and W 
must be constants and therefore a further Killing vector tangent to the Killing orbits 
occurs; for types /// and VI, W must vanish, the type /// solution then being a 
homogeneous Bianchi VI model, and the type VI one such that /i + p = 0, thus not 
corresponding to a perfect fluid. 

Notice that in []33l , the homothetic constant n was set equal to 1 from the begin- 
ning, before choosing the coordinates in the surfaces orthogonal to the Killing orbits. 
By doing so, some solutions were left out, since, although one can always re-scale the 
proper HVF X with a factor 1/n, such an scaling can not always be reabsorbed by a 
redefinition of the coordinates. We correct here that error. 

As was pointed out before, all diagonal (W = 0), perfect fluid solutions (admitting 
an orthogonally transitive Abelian G2 with flat spacelike orbits) and such that the 
metric functions A, B, R and F are separable in the variables t and z are already 



known 50, 51 



We shall next present some new exact solutions not included in ||33|| , which have 
been obtained assuming W = (diagonal), but which are not of separable variables 
in the above sense. 



Type ///: 



A 

B " 1 + 2A 



R = 1 , / = e Ht - z) , (74) 

3 (I-„)(t-*)\ 
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where < A < 1/2 and A 2 + n 2 = 1/2 and 



(/i - p)Ae n{t+z) = 



4A 



1 + 2A V2 
4A 



n 



1 + 2A 



1\ A' 



n ~2j A 



A, a and a are constants. In order to satisfy the energy-conditions, if n > 0, a must 
also be positive and the value of n is then restricted to (1/2, \[2j2)\ for negative 
values of n, a is also negative and n G (— v / 2/2, —1/2), in these cases a barotropic 
equation of state does not exist. When a = the solution is a particular case of a 
homogeneous Bianchi VI model. 



Type III: 



A(t-z) 



it> = l, / = e 

i = e (*-W-«>, A 2 = (n-1) 2 , n^i, 

B = aexp p + 2 2 n ( ^~ 1)2 (t - [exp (-(2n - l)(t - z)) - l] c 
where a and c are constants 



(75) 



» = P= 2A 



2n 



3 + c 



3 -ra(t+z) 



Type ///: 

R = 6 X ^ , / = e-^*- 2 ) , A = a 2 e (2A+2n-l)(t-r) ? £ = & 2 ? (7g) 

where a and 6 are constants, A = — (2n — 1) + y / 3n 2 — 2n + 1/4, and 
(At - p)e n( ' +2) = 4 A ( 2n - 1 - 2A) , 

£5 

(/^ + P)e^ = \^ + ±-\) + ^\{\ + 2n-\). 

From where it follows that, in order that the solution satisfies the weak and dominant 
energy conditions all over the manifold, one must have 

n E (0.8210368162407501...,^) , 

and, again, a barotropic equation of state p = p(/i) does not exist, except in the case 
n = 3/2, when /i — p. 



Type V: 



e 



t+z ,, 2c i +2c l-M^dt 2 + dz 2 \ ^ - 



18 



Po MV - 1 
^ P e'+ 2 2M 2 M 2c2 + 2c + 2 ' " ' 
where c, M and / c are constants, and ip is a function of t — z given implicitly by 

M 4 



M 2 (t-^) = \n\ip\ —ip 



2 



Type V: 



R=l, f=e 



\(t-z) 



,2(n-l)(t-«) 



5 



A 2 = ( n _ i)( n _3) , 



1 - e 



-2(n-l)(t-z) 



2(n-l) 



D - u - (^- 1 )( 4 - c ) c -na + ,) 

with a and c constants, and n restricted to being n < 1 or n > 3 in order to satisfy 
the energy conditions. 

For type VI two solutions have been obtained. For the sake of simplicity, we will 
give them in non-comoving coordinates. Thus 



ds 2 = ^ {-dt 2 + dz 2 + e~ 2t B 2 dx 2 + e- 2qt S 2 dy 2 } , 



e* e* . 

■u t = — — cosh a , m z = — smh a 
F F 



The first solution is 

-i a sin c q sin c 

F = f B, S = &-—E , 5 = C , 
H = p = e" 2t 2q(l - q)f 2 aC~ 2 " 2 ^ (A cos c + 2a sin c)E 2 , 



(79) 
(80) 

(81) 



where 



A = a 2 e^ z - e~^ z , C = a 2 e^ z + e~^ z , 



E = exp 



2q cos c 
1-9 



and 

cosh a = 



y/l + sin cC 



V 4a A cos c + 8a 2 sin c 



sinh a 



— (1 + sine) A + 2a cose 
\/l + sin cV 4a A cos c + 8a 2 sin c 



where c, a and f are constants. In order to have positive energy density, the pa- 
rameter q is restricted to q G (0, 1). A particular, simpler case can be obtained by 
choosing cose = and sine = 1. 
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The other solution is 



F = f Q B , S = R l+ ^T)T 1+ ^i , B = R^)T^ 



(82) 



-2t? 1 2 j?2 r 2 D 2 2rri2 

fi = p = e a j Q c R — q T 



o 9 n 2c r, 

R <^-<i) T~~ 



where 
and 



R = e 2 z - a 2 e * ' 



1-9 



T = e~ 2 + aV 



1-9. 



cosh a 



(c - g)Ti? 



sinh a 



gT 2 - ci? 2 



2a^/c 2 R 2 - g 2 T^ ' 2a^c 2 R 2 - q' 2 T 2 

again c, a and f are constants. Notice that the solution is only valid for 

c 2 R 2 - q 2 T 2 > . 

Notice that the solutions that have a stiff matter equation of state (p = fi) can 
be derived from vacuum solutions (also admitting an Abelian G 2 ) using a method 
proposed by Wainwright et al. |58j 



4.2 Non-Abelian G 2 

For a non- Abelian G2, a local system of coordinates can be chosen in which the Killing 

vectors, say £ and rj, are 

£ = ft , 77 = x 1 ^ + <9 2 . (83) 

Now, supposing the existence of a proper homothetic vector field, say X, one can see 
that the only allowed Bianchi type for H 3 is i.e.: 

[£,»?]=£, [£,X] = 0, fo,X] = 0. (84) 
Taking into account (|83|) and (^4|) one easily comes to the following form of X: 

X a = ^ 2 x \x\ x 4 ), X 2 (x 3 , x 4 ), X 3 (x 3 , x 4 ), X 4 (x 3 , x 4 )) . (85) 

Note that, since the orbits associated with if 3 and G2 can not coincide ||, the com- 
ponents X 3 and X 4 of the homothetic vector field cannot both vanish. Also, notice 
that the homothetic constant can in this case be set equal to one without altering 
our choice of coordinates. 



4.2.1 The orthogonally transitive case 

As in the previous Abelian case, we will restrict our attention just to the orthogonally 
transitive G 2 metrics. Regarding non-orthogonally transitive G2, a discussion can be 



found in |59[ where a study of perfect fluid solutions with four-velocity orthogonal 
to the isometric orbits is given. In that reference it is also assumed that the Killing 
vector £ is hypersurface orthogonal and the homothetic vector field, X, is orthogonal 
to the fluid velocity. These assumptions imply that the fluid is to be "stiff" (p = /i) 
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without any a priori assumption of an equation of state. No explicit solutions are 
known so far, but it is shown that solutions with pressure and matter positive on an 
open set can in principle exist by suitably specifying the initial conditions. 
For orthogonally transitive G2, the metric can be written as 



9ab 



o .2 ™2 

e an e a i2 

-X 2 

e x a i2 



V 








a-22 





\ 



^33 

ej 



where e = ±1 and = a^(a; 3 , x 4 ). 

For this case, on can see that, assuming non-null homothetic orbits V3, it is always 
possible to perform a coordinate change in the two-spaces orthogonal to the Killing 
orbits, such that it brings the homothetic vector field to the form 



X a ^( e 3;2 X 1 {x 3 ,x 4 ),X 2 {x 3 ,x 4 ),l,0) . 
and the line element can be written as 

ds 2 = A(e~ x2 dx l + Wdx 2 ) 2 + B(dx 2 ) 2 + F Udx 3 ) 2 + e(dx 4 ) 2 



17) 



where A, B, F, and W are functions of x 3 and x 4 alone. 

Specializing now the homothetic equation to the homothetic vector 
metric yields the following forms for X and the metric functions 



2 

X a = (ae x , n, 1, 0) , a, n = const 



.4 



2(l+n)x a 



a(x ) 



B 



^2x c 



b(x 4 



J2x c 



f(x 4 



and 



W 



<i_ I _ — 7! V 

n ' e 



' 3 w(x 4 ) 



—ax 3 + w(x 4 ) 



n = 



In the case n^O, one can still perform the coordinate change 



x 



,1 



a r 2 

x' e 

n 



and the 

(89) 
(90) 

(91) 
(92) 



so that the homothetic vector field and W take the forms 

X a = (0,n, 1,0) , w 



e~ nx w(x 4 ) 



(93) 



It can be easily shown, just by computing the Einstein and Riemann tensors, that 
the only vacuum solution with a non-Abelian group of isometries acting orthogonally 
transitively and admitting a proper homothetic vector field, is Minkowski space-time. 

For perfect fluid solutions, equation (|) specialized to the Killing vectors (|83|) 
implies for the fluid velocity 



u a = (e x2 u 1 (x 3 ,x 4 ),u 2 (x 3 ,x 4 ),u 3 (x 3 ,x 4 ) J U4(x 3 ,x 4 )^j 



(94) 
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For orthogonally transitive G2 models, it is possible to perform a change of coordinates 
in the x 3 , x 4 plane so as to write the four velocity and the metric as 

u a = (e~ x2 u 1 (x 3 ,x 4 ),u 2 (x 3 ,x 4 ),0,u 4 (x 3 ,x 4 )^ , (95) 

ds 2 = A(x 3 ,x 4 )(e- x2 dx 1 + W(x 3 ,x 4 )dx 2 ) 2 + B(x 3 ,x 4 )(dx 2 ) 2 

+F(x 3 ,x 4 )(dx 3 ) 2 + G(x 3 ,x 4 )(dx 4 ) 2 , (96) 

and, as a consequence, the field equations take on a much simpler form. In this 
coordinate chart and taking into account @ and (0) specified to the components gx 3 , 
9ui 923 and g 2 4, it is easy to see that X must be of the form 

X a = (ae x \ n, X 3 (x 3 ), X 4 (x 4 )) , (97) 

and the following possibilities then arise 

I X a = {ae x \n, 0,1) , (98) 

II X a = (ae x2 ,n,l,0) , (99) 

III X a = (ae x2 ,n,l,l) , (100) 

In the three cases, for one can perform the coordinate change ([92]) thus enabling 
one to set a zero. 

The homothetic equation specialized to the metric (|96| ) yields then the following 
possibilities 



type 




A(x 3 , x 4 ) 


W(x 3 ,x 4 ) 


La 


(0,n,0,l) 


e 2{l+n ^ a{x 3 ) 


e _nx4 w(x 3 ) 


1.6 


(ae x \ 0,0,1) 


e 2x4 a(x 3 ) 


—ax 4 + w(x 3 ) 


Il.a 


(0,n,l,0) 


e 2{1+n ^a(x 4 ) 


e~ nxi w(x 4 ) 


II.6 


{ae x \ 0,1,0) 


e 2x6 a{x 4 ) 


—ax 3 + w(x 4 ) 


IILa 


(0,71,1,1) 


e (l+n)(x a +x' l ) a / x 3 _ x 4\ 


e~^ < " xi+x4 '^w{x 3 — x 4 ) 


III.6 


(ae x \ 0,1,1) 


e x6+x *a(x 3 — x 4 ) 


-f(x 3 + x 4 ) +w(x 3 -x 4 ) 



Table 4.5 



type 


B(x 3 , x 4 ) 


F(x 3 , x 4 ) 


G(x 3 ,x 4 ) 


I 


e 2xll b(x 3 ) 


e 2xi f(x 3 ) 


ee 2x * f(x 3 ) 


II 


e 2x %(x 4 ) 


e 2xa f(x 4 ) 


ee 2xA f(x 4 ) 


III 


e xJ+x ' t b(x 3 — x 4 ) 


e xi+x ' l f(x 3 — x 4 ) 


e xi+xi g(x 3 — x 4 ) 



Table 4.6 
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and for the velocity field one has 



type 


Ui{X , X j 


„~, ™4\ 

■u 2 (^ °,X ) 


U4(X , X J 


T n 


Jl+n)x* f /„3\ 


e 1*2 v-c j 


^3 ^ 01 . ( n™ 3 ^ 


Lb 


e xi Ui(x 3 ) 


— aa^e* -u^a; 3 ) + e x ii2(x 3 ) 


e^iL^a; 3 ) 


ll.a 


e^+^u^x 4 ) 


e x *ii2(x 4 ) 




ll.b 


e xi u\{x A ) 


—ax 3 e x u\(x 4 ) + e x ti2(x 4 ) 




III.o 


1+ " (x 3 1 X 4 ) " 1 \ 

e 2 \ x +x >u\[y) 


e 2 u 2 (v) 


e 2 -u 4 (t>) 


111.6 


e 2 -u^wj 


—a ^ e 2 -uxfvj + e 2 u 2 (vj 


X^ -\~X^ A / \ 

e 2 m 4 (wJ 



Table 4.7 



where the subcases a and 6 refer to whether n 7^ orn = respectively, and v stands 
for x 3 — x 4 . 

Since 113 = in this coordinate chart, the components of the Einstein tensor G13, 
G23 and G34 must vanish identically, hence 



W* = , 







B. 



IA4 1^4 

U 2 F 



A3 

3 5 4 
4~1T 



^3 

1^3 /Ai 
4 G V A 



B 



B 
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(101) 
(102) 
(103) 



BV4A IF 4 B J 4 G \ A B J B 

where ^ means derivative with respect x l . From these equations we found more ex- 
plicit forms for the metric functions, that are given by: 



Case I. a 

ds 



2 2x4 > k e 2nxA 



f(x 3 
2 



j-j-rt 2 2-j-n 

" (e^^+we-" 14 ^ 2 ) + &[/(x 3 )] n (da; 2 ) 



+f(x 3 ) ( (dx 3 ) + e (dx 4 ) 



(104) 



where /c = ±1, w and 6 are arbitrary constants. Note that n must be different from 
zero or a third Killing vector occurs and the metric becomes then LRS. In such a 
case the metric functions a(x 3 ) and b(x 3 ) are still proportional to each other, but no 
relation exits, in principle, with f(x 3 ). 



Case I.b 



ds 2 = e 2x4 ja(x 3 ) (e'^dx 1 - ax 4 dx 2 ^ + ba(x 3 ) (dx 2 ) 2 + (dx 3 ) 2 + e (dx 4 ) 2 J , 

(105) 

where b is an arbitrary constant different from zero to prevent a non-singular metric, 
and a 7^ if the group G 2 is to be maximal. 
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Case II 

Only one possibility arises in this case, namely a = n = 0; thus the line element 
becomes 

ds 2 = e 2x3 L(x 4 ) (e~ x2 dx l + w(x 4 ) dx 2 f + b(x 4 ) (dx 2 ) 2 + (dx 3 ) 2 + e (dx 4 )'' 



(106) 
Case III. a 

ds 2 = e x3+x4 i [ e n ^ + ^a(v)(e- x2 dx 1 + we- nx4 dx 2 ) 2 + 

be n{x3 - x4) a(v) (dx 2 ) 2 + f{v) (dx 3 ) 2 + g{v) (dx 4 )^ , (107) 

w being an arbitrary constant and v = x 3 — x 4 . 
Case IH.b 

\e- x2 dx 1 - ax 4 dx 2 f + b [dx 2 ) 2 + f(v) (dx 3 ) 2 + g{v) (dx 4 )' 



ds 2 = e x3+x4 1 a(v 



(108) 

Again, a and n are arbitrary non-null constants. The differential equation ( |103| ) for 
case III can be rewritten as 

o^f^wf^^u^^y+^v, (109) 

9 V a I J \ a J a \a J \a J 

where the prime denotes an ordinary derivative with respect to the variable v = 

3 4 
X — X . 

4.2.2 Diagonal case 

Since the field equations for a perfect fluid are still so complicated we will make a 
further assumption that will bring the metric into diagonal form; namely: the Killing 
vector £ being hypersurface orthogonal. 

The possibilities are now restricted just to diagonal subcases (a), since diagonal 
subcases (b) do always admit a further Killing vector tangent to the Killing orbits V2. 

Computing the Einstein tensor for those metrics, one has 

G u = and G 24 ± . (110) 

Consequently, we have chosen a coordinate chart in such a way that the fluid flow 
velocity always lies in the two plane spanned by d/dx 2 and d/dx 4 at each point. 
Therefore, we will have in all cases 

u = u 2 dx 2 + u,dx 4 , ^£ + ^ = -1 . (Ill) 

922 #44 
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A careful analysis of all the possibilities reveals that, in most cases, there exist 
further KVs and, in the instance of null homothetic orbits, the energy conditions 
cannot be fulfilled. Apart from these cases, it is worth mentioning that whenever X 
is orthogonal to u, the metric and field equations are 

ds 2 = e 2x " jaVK 2 * 2 (dx 1 ) 2 - eb 2 (x 4 ) (dx 2 f + (dx 3 f + e (cfe 4 ) 2 } , (112) 








1 a'b' a" 
(r ab a 



'a' 


b'~ 


2 


'b" 


a"' 




a 


b_ 




b 


a 





1 b l 

'b 2 + ab b 



(113) 
(114) 



a prime indicating a derivative with respect to (x 4 ), and then one necessarily has 



p = H = e 



-2x A 



-4' 



(115) 



without previously assuming the existence of a barotropic equation of state. 
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